SHEP 93/94-19 
hep-ph/9406303 

Unification Constraints in the Next-to-Minimal Supersymmetric 

Standard Model 



^ ■ T. Elliott, S.F. King * 

^ ' Physics Department, University of Southampton, 

Southampton S09 5NH, UK. 

6 

\D : P.L. Whitet 



Theoretical Physics, University of Oxford, 



in 

> 

O ■ 1 Keble Road, Oxford 0X1 3NP, UK. 

m 

^ '. (February 1, 2008) 



Abstract 



^ . We perform a renormalisation group analysis of the next-to-minimal su- 

persymmetric standard model (NMSSM) based on the following constraints: 
two- loop gauge coupling unification at a variable scale Mx, running the 
gauge couplings through the low energy thresholds; universal soft supersym- 
metry breaking parameters; correct electroweak symmetry breaking. The 
phenomenological implications of our results include a standard model-like 
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Higgs boson of mass in the range 70-140 GeV. 



Typeset using REVT^ 



Supersymmetry (SUSY) [|l| is an attractive candidate for new physics beyond the stan- 
dard model at the TeV scale. The existence of SUSY would greatly ease the naturalness 
problems associated with the construction of a grand unified theory (GUT) of the strong 
and electroweak interactions. The basic idea of GUTs is that the gauge couplings, which 
govern the strength of the strong and electroweak interactions at low energy, are actually 
equal to some unified coupling gx at some very high scale Mx due to their renormalisation 
group (RG) running 0. The original motivation for SUSY broken at a TeV was to help 
to stabilise the Higgs mass against GUT scale quadratic radiative corrections. The combi- 
nation of SUSY and GUTs has recently found some indirect experimental support due to 
the accurate measurement of the strong and electroweak couplings on the Z pole by LEP. 
These measurements are inconsistent with gauge coupling unification if a standard model 
desert is assumed, but are consistent with unification if a SUSY desert above the TeV scale 
is assumed 0. 

In the first post-LEP analyses the whole SUSY spectrum was either assumed to 
be degenerate at the scale Msusy, or smeared around this scale. In reality the SUSY 
partners may have a complicated spectrum, parametrised by a large number of soft SUSY- 
breaking parameters, which may spread over one or two orders of magnitude of masses. In 
order to reduce the number of independent soft SUSY-breaking parameters one may appeal 
to supergravity or superstring scenarios which involve the notion of SUSY breaking in a 
hidden sector coupled only gravitationally to our observable sector. By this means, in the 
minimal supersymmetric standard model (MSSM), one ends up with just four independent 
soft SUSY breaking parameters: mo, M1/2, Aq and Bq corresponding to the universal soft 
scalar, gaugino, trilinear and bilinear couplings, respectively. Using these four parameters, 
together with the top quark Yukawa coupling which plays an important role in driving 
electroweak symmetry breaking several groups have performed an RG analysis whose 
goal is to predict the SUSY and Higgs spectrum and then use this spectrum as the basis 
of a more reliable estimate of gauge coupling unification by running the gauge couplings 
through the various thresholds from the Z mass Mz up to a unification scale Mx [Hll- 
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Apart from the constraints of gauge coupling unification and correct electroweak symmetry 
breaking, there are various other phenomenological and cosmological constraints which may 
be apphed and recent studies have concluded that it is possible to satisfy all these constraints 
simultaneously 0. 

Here we shall consider a slightly different low energy SUSY model, but one which is 
equally consistent with gauge coupling unification, namely the so-called next-to-minimal 
supersymmetric standard model (NMSSM) 10-^. The basic idea of the NMSSM is to add 
just one extra gauge singlet superfield to the spectrum of the MSSM, and to replace the 
/i-term in the MSSM superpotential with a purely cubic superpotential, 

fxHiH2 \NH^H2 - (1) 

The motivation for this "minimal non-minimal" model is that it solves the so-called /i- 
problem of the MSSM in the most direct way possible by eliminating the /z-term alto- 



gether, replacing its effect by the vacuum expectation value (VEV) < >= a;, which may 
be naturally related to the usual Higgs VEVs < Hi >= Ui. However there are other solutions 
to the yU-problem |TD|. Also, the inclusion of singlets may cause the destabilisation of the 



hierarchy if there are strong couplings to super-heavy particles such as Higgs colour triplets 
|Tl[. Recently similar effects have been shown to result from non-renormalisable operators 



suppressed by powers of the Planck mass . The dangerous non-renormalisable operators 
would require that gravity violate the symmetry which is respected by the renormalis- 
able operators of the theory. Our view is that since these effects are model-dependent the 
NMSSM is as well motivated as the MSSM and should be studied to the same level of ap- 
proximation. Only by so doing may the two models be phenomenologically compared when 
(or if) Higgs bosons and SUSY particles are discovered. 

In this Letter we present our first results of a GUT-scale RG analysis of the NMSSM, 
imposing simultaneously the constraints of correct electroweak symmetry breaking and cou- 
pling constant unification. In our analysis we shall drop all quark and lepton Yukawa 
couplings apart from that of the top quark ht- We shall assume a soft SUSY breaking 



potential which gives rise to three universal soft parameters at the Mx, namely mo, M1/2, 
Aq, mentioned above 0. We shall use the two-loop gauge and Yukawa SUSY RG equations 
and one-loop soft SUSY RG equations to integrate from the GUT scale Mx to a low energy 
renormalisation scale Q. Below Mx, the three soft mass parameters evolve into 32 separate 
soft parameters, corresponding to 3 gaugino masses Mj, 11 trilinear couplings Ai (those of 
the MSSM plus Ax and Ak) and 18 scalar masses (including the soft Higgs masses m^., 
m^). We assume no inter-generational mixing, and since we drop Yukawa couplings which 
are small compared to ht, the first and second generation soft masses will run identically. 
The unification constraints are: gf{Mx) = g'x, X{Mx) = Aq, k{Mx) = ko, ht{Mx) = hto, 
Mi{Mx) = Ml/2, m^{Mx) = ml Ai{Mx) = Aq. 

In practice we shall scale all dimensionful quantities by appropriate powers of M1/2, and 



denote the resulting dimensionless quantities by a tilde. At the low energy MS scale Q we 
find the global minimum of the one- loop effective neutral Higgs potential, scaled by M1/2, 
including loops of top quarks and stop squarks [0,|l3 • We only consider the neutral scalar 
Higgs potential, since there are simple analytical conditions which test for the presence of 
squark and slepton VEVs Although these are not guaranteed to prevent all possible 
charge or colour breaking minima as discussed by Gunion et al in reference they are 
sufficient to avoid the most dangerous case of slepton VEVs, while a complete analysis of 
all possible colour breaking minima in the one loop effective potential is not feasible. We 
implement these conditions after running to an energy scale of order a typical slepton VEV. 
In addition we ensure that the mass squared of the charged Higgs remains positive in order 
to avoid breaking electromagnetism. The result of minimisation yields the scaled VEVs 
Ui, X, from which the condition u = 174 GeV fixes Mi/2- Given M1/2 this enables all the soft 
parameters and the VEVs to be determined and, as described elsewhere p!3|JI^ , it is then 



straightforward to calculate the top and stop corrections to the Higgs boson mass spectrum. 



^ Note that the B parameter does not occur in the NMSSM since we have set ^ = 0. 
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The SUSY spectrum is calculated at tree- level using standard results [|T[, apart from 
squark masses which we calculate using the one-loop effective potential. There are a number 
of phenomenological constraints which will cut down the allowed regions of parameter space. 
We shall require that all sleptons and stops are heavier that 43 GeV and all charginos 
are heavier than 47 GeV. Charged Higgs bosons are required to be heavier than 45 GeV. 
Gluinos and squarks other than stops are required to be heavier than 100 GeV. The lightest 
neutralino is required to be the lightest SUSY particle. The lightest CP-even neutral Higgs 
boson h is required to satisfy mh/Rzzh ^ 60 GeV, where is the mass and Rzzh is the 
ZZh coupling scaled by the standard model coupling. We shall require a large top quark 
mass, rrit > 150 GeV [1^, which puts a strong restriction on the parameter space of this 
model. In our analysis we input the gauge couplings gi{Mz) and g2{Mz) and run them up 
through the 24 SUSY and Higgs thresholds to find Mx and gx, which must of course be 
consistent with our original input values. By iterating this procedure we obtain solutions 
which satisfy both the requirements of correct electroweak symmetry breaking and coupling 
constant unification simultaneously. According to our procedure the parameters Mx and 
gx are determined from gi{Mz) and g2{Mz), and so our input parameter set is: Aq, kQ, hto, 
TTiQ, Aq, with the value of M1/2 determined from the VEV. The physically relevant quantities 
as{Mz) = gz^{Mz) I^Ti, mt, r = x/u and tan/? = 1^2/1^1 are all outputs 0. 

A somewhat similar analysis to that described has been performed in the NMSSM by 



another group [lf:-18|. We have properly included the low energy threshold effects in our 



^Our quoted values of mt always refer to the one-loop physical pole mass mt = m^° ^ = mt{mt)\\ + 
^as(mi)]. Also, a word about our conventions is in order. In the convention of Eq.(l), we input 
positive values of Aq, A;o, and both positive and negative values of Aq, and look for solutions with 
both positive and negative VEVs, except that we require V2 > 0. Such solutions can always be 
re-interpreted in terms of purely positive VEVs using the three symmetries of the one-loop effective 
potential: (i) A —X,k —k,x — x, (ii) A — A,z^i —vi, (iii) J^i — > —1^1,1^2 —1^2- 
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unification analysis, whereas in ref. [|T6|-0 such effects were ignored. Since our choice of 
0:3 (M^) and of Mx and gx are determined by the requirements of consistent unification and 
thus vary with differing spectra, our approach is consequently much more computationally 
intensive since we have to run down the couplings, minimise the potential, calculate the 
spectrum, then run the couplings up through the thresholds, iterating this procedure many 
times in order to obtain a single point in parameter space which is consistent with unification 
and all the other constraints. 

We initially considered a grid of 360,000 values of parameters in the ranges: Aq = 0.01 — 
2.0, = 0.01-2.0, hto = 0.5, 1,2,3, rho = 0.2-5.0, Ao/mo =-4 to +4, then later investigated 
the phenomenologically interesting regions in more detail. For the acceptable regions of 
parameter space we find that rrit > 150 GeV is attainable for ko < Xo < hto, with Aq < 1 
and 0.5 ~ hto ~ 3, where the upper limit on hto is a requirement of perturbation theory. 
Although all values of hto in this range lead to consistent solutions with | Ao/mo | ~ 3 — 4, 
for hto ~ 0.5 we also find solutions for smaller lAo/mol, down to zero. However the cut 
on the top mass of rrit > 150 GeV excludes this region, and so we discuss it no further 
here. The successful values of mo vary across the whole range, but are correlated with hto 
and Ao/mo- To be specific, we find for |Ao/mo| ~ 3 that hto = 0.5, 1,3, is associated with 
mo ~ 0.5 — 1, 1 — 2,2 — 5, while smaller values of mo can be obtained for larger |Ao/mo|, 
although the maximum value of | Ao/mo | is about 4 (because of the slepton VEVs constraint). 
There are also correlations between Ao and ko at the unification scale, as discussed later. 
In the successful regions of parameter space we find | tan/?| ~ 3 — 20 (or larger, but our 
approximations break down here), while |r| is virtually directly proportional to M1/2, with 
|r| 3> 1 even for small Mi/2- Thus, as in ref. |T^, we find that the allowed region of parameter 
space is characterised by small dimensionless couplings |Ao| < 1, |A;o| < 1 and large singlet 
VEV 1^1 > 1000 GeV, corresponding to the approximate MSSM limit of the model 0. 

To illustrate the correlation between Ao and ko, in Fig.l we show contours of hto in the 
Ao — ko plane with mo = 2, M1/2 = 500 GeV and |y4o/mo| = 3. Although we could have 
chosen other values of M1/2, the qualitative features of Fig.l would be unaltered, and the 
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value of ko would only change by a few per cent, since in this region the value of M1/2 is a very 
sensitive function of Aq and k^. The htQ = 1, 2, 3 contours correspond to rrit ~ 175, 185, 190 
GeV, respectively, with | tan/3| varying from 3-7 and |r| from 100-30, from the lowest to the 
highest value of ko respectively. Values of ko beyond the ends of the plot are forbidden by the 
requirements of electroweak symmetry breaking, although the exact ranges of acceptable ko 
could be altered by varying our choice of Aq and mo slightly. Similarly, nit would be altered 
if we use a different M1/2 primarily due to the resulting change in as{Mz). As observed by 
other authors we find that larger values of M1/2 result in smaller values of as{Mz). It 
should be noted that although changing the sign of Aq leads to roughly the same output 
data with a change in the sign of r, this is only approximate, and in particular one sign of 
Ao leads to a less restrictive slepton VEV constraint. 

We find that, for many choices of our input parameters, there appears to be a certain 
minimum value of M1/2 (recall that this is an output parameter). For example, if we impose 
a cut of Ml/2 ~ 150 GeV, then rrit > 150 GeV is only possible for hto ~ 1.5 — 3, |y4o/r?7,o| ~ 3, 
rho 3 — 5, corresponding to rrit ~ 185 — 195 GeV. The origin of the minimum value of M1/2 
is due to a different competing global minimum of the one-loop effective potential, which 
becomes preferred below a minimum M1/2 value. This happens because the logarithm in 
the one-loop effective potential depends both on the renormalisation scale Q (which we have 
so far taken to be 150 GeV) and on Mi/2- To test the reliability of our results, we have 
investigated their Q-dependence. In doing this we have consistently used the same value of 
Q both in the RG running of masses and couplings and in the one-loop effective potential, 
since there is a significant cancellation between the effects of changing Q in these two places. 

Although varying Q leads to very little change in the VEVs for any given minimum, it 
can alter the VEVs substantially by changing which of the minima is deepest (there are 
often minima with very different VEVs but with similar values of the potential). Given 
this behaviour, one must be very careful about drawing overly restrictive conclusions based 
on which of the various minima is preferred for a given region of parameter space. We 
find that reducing Q will typically lead to a reduction of the minimum value of M1/2 in 
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those regions where small M1/2 is not possible for Q = 150 GeV, but does not otherwise 
qualitatively change the behaviour. Note that the range of parameters htQ ^ 1.5 — 3, 
l^o/'^ol ~ 3, rfiQ ^3 — 5 appears to be a "safe" range from the point of view of achieving 
small values of M1/2 ~ 150GeV independently of the choice of Q. This dependence on Q 
is an inherent problem in this model because of the form of the effective potential which 
has several competing minima. In both our analysis and that of ref. [[T6| -|T^ electroweak 
symmetry breaking has been analysed using the one-loop effective potential, including top 
quarks and stop squarks. In principle the Q dependence could be eliminated by including 
all particles in the effective potential, and possibly working beyond one-loop. Here we shall 
restrict ourselves to a discussion of the Q-dependence of our results. 

We now consider two examples of the sort of Higgs and SUSY spectrum one might 
typically expect in the constrained NMSSM, the first with relatively large values of hto and 
mo, and the second with relatively small values of hto and rfiQ. These two extreme cases 
span the complete range of spectra predicted by the constrained NMSSM, subject to the 
constraint of a large top mass > 150 GeV. 

In Fig. 2 we plot the spectrum as a function of M1/2 for the case hfo = 2, |Ao/mo| = 3, 
mo = 5, Aq = 0.4. With Q = 150 GeV we find that M1/2 niay go down to small values, which 
is not unexpected since these parameters are in the "safe" range hto ~ 1-5 — 3, |Ao/mo| ~ 3, 
rho ~ 3 — 5, referred to above. Note that the experimental constraint that the charginos are 
heavier than 47 GeV implies that M1/2 ~ 70 GeV in this case. For M1/2 ^ 100(1000) GeV, 
which is controlled by choosing ko = 0.275(0.300), we find |tan/3| ~ 6(8), |r| ^ 13(120), 
as{Mz) ^ 0.121(112), Mx = 2.1(0.99) x 10^^ GeV, gx = 0.71(0.69). The lightest CP-even 
Higgs boson has standard model-like couplings {Rzzh > 0.99 everywhere) and for M1/2 ~ 100 
GeV is in the LEP2 range, as are the lighter chargino and neutralinos which have a much 
stronger M1/2 mass dependence. The top quark mass ranges from m^ = 193 — 184 GeV, 
being smaller for larger M1/2 due to as{Mz) being consequently smaller. For M1/2 ~ 100 
GeV, the lightest stop and gluino are not too much heavier than the top quark, although 
the remaining sparticles and Higgs bosons are significantly heavier than the top quark. 
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In Fig. 3 we plot the spectrum as a function of M1/2 for the case hto = 0.5, |Ao/mo| = 
4, mo = 0.5, Ao = 0.1. These parameters are outside the "safe" range hto ^ 1.5 — 3, 
l^o/'^ol ~ 3, rho ^ 3 — 5, in which small M1/2 can be achieved independently of the choice 
of the renormalisation scale Q, and so we plot the spectrum for two choices of Q. The 
main effect of changing Q is to change the value of M1/2 at which the data cuts out. For 
Q = 150(25) GeV, we find the minimum values M1/2 = 300(125) GeV corresponding to 
mo = 150(62.5) GeV. As in Fig. 2, the lightest CP-even Higgs boson has standard model-like 
couplings, and for a given M1/2 is even lighter than in Fig.2. For the smaller M1/2 values 
which we can obtain with Q = 25 GeV the lighter chargino and neutralinos may be in the 
LEP2 range. Unlike Fig.2, the left-handed sleptons are now much lighter (due to the smaller 
value mo = 0.5) while the lighter stop is much heavier (due to the smaller value hto = 0.5). 
The top quark mass has a maximum value of m^ = 175 GeV for M1/2 = 125 GeV and here 
I tan/?| = 3.4. The gluino in Fig.3 is now the heaviest sparticle, whereas in Fig.2 it was one 
of the lighter ones. In general the Higgs and sparticle masses in Fig.3 are focussed into a 
narrower band of masses than in Fig.2, which is a simple result of having mo = 0.5 rather 
than mo = 5. 

Finally it is worth comparing our spectra to the spectra discussed in a paper which 



appeared just as this present article was being finalised [^. In this most recent analysis 
the Higgs and SUSY spectrum was presented as as series of scatter plots for masses and 
couplings. Most of this data and all of ours corresponds to one of the CP-even Higgs bosons 
being almost pure gauge singlet. Such a decoupled Higgs boson might be expected since the 
constrained NMSSM is close to the MSSM limit. In the analysis of ref. [|18| the decoupled 



Higgs boson may have a mass either less than or greater than the lighter of the other two CP- 
even Higgs bosons and when the would-be decoupled Higgs is close in mass to the lighter 
physical Higgs, strong mixing can occur, leading to two weakly coupled Higgs bosons of 
the kind discussed in ref. |T^. However, given the range of parameters considered in this 
paper, we find that the decoupled CP-even Higgs boson is always substantially heavier than 
the lighter of the two physical CP-even Higgs bosons. Similarly the CP-odd Higgs bosons 

10 



are much heavier than the hghtest CP-even Higgs boson. The reason for this difference is 
simply that in ref. the range of parameters considered exceeds the range considered 
in this paper. It turns out ^ that in order to bring down the mass of the CP-even singlet 
sufficiently one requires rho <^ 0.2, Aq < 0.01 and ko <^ Aq, corresponding to extremely large 
values of r ^ 100. We have checked that, in this region of parameter space, the singlet 
CP-even Higgs boson does indeed become much lighter, leading to the strong mixing effect 
mentioned above. However, for the range of parameters considered in the present paper, a 
standard model-like Higgs boson with a mass in the range 70-140 GeV is preferred. 
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Figure Captions 



Figure 1: A contour plot showing allowed values of Aq and ko for given values of 
hto = 1, 2, 3 (going from left to right) in the Aq-Zcq plane, corresponding to M1/2 = 500 GeV, 
rho — 2 with |Ao/mo| = 3. 

Figure 2: Masses of particles as a function of M1/2. The input parameters are hto — 2, 
\Ao/mo\ = 3, mo = 5, Aq = 0.4, ko = 0.275 - 0.3. Wc use Q = 150 GeV. Neutralinos (solid 
lines) , charginos (dot-dashed lines) , CP-even Higgs (short dashed lines) , lighter stop and top 
quark (both dotted hues), left-handed sleptons (long dashed lines), and gluino (quadruple 
dashed hnes) are displayed. One of the CP-odd Higgs and the charged Higgs bosons (not 
shown) are roughly degenerate with heavier CP-even Higgs bosons. The remaining CP-odd 
Higgs boson (not shown) is heavier than the heaviest neutralino. The lightest CP-even Higgs 
(which has standard model- like couplings), lighter chargino and lightest two neutralinos are 
all in the LEP2 range for M1/2 ^ 100 GeV. 

Figure 3: Masses of particles as a function of M1/2. The input parameters are hto — 0.5, 
\Ao/mo\ = 4, mo = 0.5, Aq = 0.1, ko = 0.07 - 0.13. In addition to using Q = 150 GeV, 
as in Fig. 2, we also show results for Q = 25 GeV, corresponding to minimum values of 
Ml/2 — 300 GeV and M1/2 = 125 GeV, respectively. Neutralinos (solid hnes), charginos 
(dot-dashed lines), CP-even Higgs (short dashed lines), lighter stop and top quark (both 
dotted lines), left-handed sleptons (long dashed lines), and gluino (quadruple dashed lines) 
are displayed. One of the CP-odd Higgs (not shown) is roughly degenerate with the second 
heaviest CP-even Higgs boson. The remaining CP-odd Higgs boson (not shown) is roughly 
equal to the lighter stop mass. The lightest CP-even Higgs (which has standard model-like 
couphngs as in Fig. 2) has a mass ~ 100 GeV. 
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